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Abstract 

A particularly simple Lorentz-violating modification of the Maxwell theory of photons maintains 
gauge invariance, CPT, and renormalization. This modified-Maxwell theory, coupled to standard 
Dirac particles, involves nineteen dimensionless "deformation parameters." Ten of these parameters 
lead to birefringence and are already tightly constrained by astrophysics. New bounds on the 
remaining nine nonbirefringent parameters can be obtained from the absence of vacuum Cherenkov 
radiation for ultrahigh-energy cosmic rays (UHECRs). Using selected UHECR events recorded at 
the Pierre Auger Observatory and assigning pseudo-random directions (i.e., assuming large-scale 
isotropy), Cherenkov bounds are found at the 10~^^ level, which improve considerably upon current 
laboratory bounds. Future UHECR observations may reduce these Cherenkov bounds to the 10^^^ 
level. 
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I. INTRODUCTION 



It is possible that entirely new phenomena at the high-energy frontier of elementary 
particle physics (with characteristic energy Epianck = ^hcyC ^ 1.2 x lO^^ GeV) lead to 
Lorentz-violating effects in the low-energy theory, in particular, the Maxwell theory of pho- 
tons. The simplest CPT-invariant Lorentz-violating modification of the quadratic Maxwell 
action involves nineteen dimensionless parameters [l|, Q, 0] . 

Existing bounds on these nineteen "deformation parameters" of modified- Maxwell the- 
ory are as follows: ten birefringent parameters are bounded at the 10~^^ level or better from 
astrophysics (see Refs. 0, [sl and references therein) and nine nonbirefringent parameters 
are bounded at the 10~^^ level or worse from laboratory experiments (see Refs. p, 0, 0] and 
references therein). 

Following earlier suggestions j^, [lo|, it has been noted [ll|, [l^ that ultrahigh-energy 
cosmic rays (UHECRs) have the potential to place further limits on these nonbirefringent 
parameters by the inferred absence of "vacuum Cherenkov radiation." In fact, Ref. [l^ 
already gave a bound at the 10~^^ level on one nonbirefringent parameter from the existing 
UHECR data (cf. Refs. [isl, Q). In this article, we obtain even better bounds for all 
nine nonbiref rinp: ent parameters from selected UHECR events recorded by the Pierre Auger 
Observatory [15] over the period January 2004 to February 2006 (l6| . 

The present article is organized as follows. In Sec. [Tll we review the theory considered, 
quantum electrodynamics with a modified-Maxwell term, and establish our notation. Also, 
we recall the Cherenkov threshold condition on which our bounds will be based. In Sec. Illlt 
we discuss the selection of an appropriate set of UHECR events from two years running 
of the partially completed Auger experiment. In Sec. llVt we obtain from the Cherenkov 
threshold condition and the selected Auger events (with pseudo-random directions) the 
promised bounds on the nine deformation parameters of nonbirefringent modified-Maxwell 
theory. In Sec. |Vl we summarize our findings and compare these new astrophysics bounds 
with previous laboratory bounds. 

Throughout this article, we employ Cartesian coordinates [x'^) = (x°, x) = {ct, x^, x^, x^) 
and the Minkowski metric {rj^^) = diag(+l,— 1,— 1,— 1). Indices are lowered with the 
Minkowski metric ?7^jy and raised with the inverse metric rj'^'^. Repeated upper and lower 
indices are summed over (Einstein summation convention). The direction of a 3-vector x is 
given by the unit 3-vector x = x/\x\. Natural units with c = h = 1 are used, unless stated 
otherwise. 
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II. THEORY 



A. Nonbirefringent modified-Maxwell theory 

The action of modified-Maxwell theory 0, [sl is given by 



modM 



1 (^/.P^ 



ua _|_ j^fivpa 



(1) 



where F^^{x) 



df,A^{x) 



diyA^{x) is the standard Maxwell field strength of gauge fields 
An{x) propagating over fiat Minkowski spacetime with a metric r/^j/ as defined in Sec. HI Here, 



corresponds to a constant background tensor (with real and dimensionless components) 
having the same symmetries as the Riemann curvature tensor and a double trace condition 
K^^p^ = 0, so that there are 20 — 1 = 19 independent components. All components of 
the K-tensor in ([1]) are assumed to be very small, <C 1, in order to ensure energy 

positivity The photonic action term ([T]) is gauge- invariant, CPT-even, and power- 
counting renormalizable (cf. Ref. (l?!). 



In order to restrict modified-Maxwel 
Ansatz for k^'^p'^ has been suggested [s]. 



theory to the nonbirefringent sector, the following 
l3: 



1 ^pp ~va _ ^vp 



in terms of the nine components of a symmetric and traceless matrix n 



pLV 



K 



up 



0. 



(2) 



(3) 



The nine Lorentz- violating "deformation parameters" k^" , called "coupling constants" in 
Ref. 



12| . can be rewritten as follows: 



{k^^") = diag(l, i, |, I) + {6k'''') , (5^00 = , 



(4) 



with one independent parameter kP^ for the spatially isotropic part of k'^" and eight inde- 
pendent parameters in SH'^'^. For later use, we already define 
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where a denotes the corresponding vector in the parameter space with squared Euchdean 
norm 

(6) 

1=0 

At this moment, it can also be mentioned that, for modified- Maxwell theory, the phase 
velocity is given by Vphasc = ka;/|k| = ck(l — S(k)), with a dimensionless function S 
depending on the wave vector dire;tio„ k aud the compoueuts of the .-tensor 0. Th. 
result implies that, for modified-Maxwell theory, the phase velocity Vphase and the front 
velocity Vfront = lini|k|^oo Vphase(k) are equal for each direction k. Moreover, the phase and 
group velocities are equal to leading order in the Lorentz-violating parameters, Vgroup = 
duj/d\<i = Vphasc + O(k^). Recall that the front velocity is the relevant quantity for signal 
propagation and causahty and that the group velocity at the dominant frequency component 



of a broad wave packet gives the velocity of energy transport; see, e.g., Refs. [l8|, [l9|. 

We now add spin-i particles with electric charge e and mass M to the theory, taking 



the usual minimal coupling to the gauge field [l7j. The action of these charged particles is 
assumed to be Lorentz invariant, so that, for the particular spacetime coordinates employed, 
the Lorentz violation of the combined theory resides solely in the k term of ([T]) . Specifically, 
the action with the standard Dirac term included is given by 

5modM+standD = ^modM + / d'^X Ipix) ('^^'{i 8^ - eA^{x)) - M)iIj{x) . (7) 

Note that the maximal attainable velocity of this charged particle equals c, which may 
or may not exceed the phase velocity of light discussed in the previous paragraph. See 
Refs. 



20 



21 



22| for further discussion on microcausality in Lorentz-violating theories and, 
in particular, Eq. (50) of Ref. [21] for microcausality in the fermionic coo^model, which is 
formally related to a special case of modified-Maxwell theory ([7]) by an appropriate linear 
coordinate transformation 0.0. 

As mentioned in the Introduction, new phenomena at the energy scale ii^pianck ~ 10^^ GeV 
may lead to Lorentz violation in the low-energy theory, possibly described by modified-Max- 
well theory ([1]). The crucial point to realize is that, a 'priori^ the Lorentz-violating parameters 
^iivpa ^gg(j j^Q^ small (e.g., suppressed by inverse powers of -Epianck) but can be of order 



unity |23|, l2J, |25| . For this reason, it is of importance to obtain as strong bounds as possible 
on all deformation parameters k^^'"^ . 



B. Cherenkov threshold condition and UHECR bounds 



Vacuum Cherenkov radiation for modified quantum electrodynamics ([7]) wit 
tonic action ([I])-([2]) has been studied in the classical approximation by Altschul 



1 the pho- 
11| and at 



4 



tree level by Kaufhold and Klinkhamer [12|. As explained in Ref. [l2j (see also Ref. |10|). 
the radiated energy rate of a primary particle with point charge Z^j-i^e, mass Mprim > 0, 
momentum ^pj-im, and ultrarelativistic energy -Eprim ~ c l^priml is asymptotically given by 



dW^modM(9prim) ^2 N /», /0\ 
~ ^prim ^ '5 l^prim j %nm/ " ' 1° j 

where the nonnegative dimensionless coefficient ^ results from appropriate contractions of 
the K-tensor ([21) with two rescaled g-vectors. The coefficient ^ may or may not depend on 
the primary particle direction qpi-jm- 

The asymptotic behavior shown in ([8]) holds only for particle energies -Eprim well above 
the Cherenkov threshold, which has the following order of magnitude: 



-E^thresh ~ AfpHm C / V K , (9) 

for an appropriate scale k obtained from the H^^'^ components written in terms of the pa- 
rameters a' (the scale k is effectively set to zero if Cherenkov radiation is not allowed). For 
completeness, we mention that, at present, the most detailed study of vacuum Cherenkov 
radiation has been performed for another Lorentz- violating theory, Maxwell-Chern-Simons 

1 for further discussion and references. 



theory, and refer the reader to Refs. [12, l26l. |27, 



28 



Continuing with nonbirefringent modified-Maxwell theory, it is now possible to derive 
an upper bound on a particular combination of the deformation parameters a' from the 
observation of a single UHECR event with a nuclear primary moving in the direction qprim 
and having an ultrarelativistic energy Ep_ » Mp_ The argument is remarkably simple 
[9|, [10| : an UHECR proton or nucleus can arrive on Earth only if it does not lose energy by 



vacuum Cherenkov radiation and this requires the particle energy -Eprim to be at or below 
threshold, Ep,-,^, < ^^thresh- 

The caveat of this simple argument is that the radiation rate ([8]) should not be suppressed 
by an extremely small numerical factor entering the coefficient ^. But there is no reason 
to expect the presence of such an extremely small numerical factor; cf . Refs. (lol . [l^ . In 
fact, a recent calculation 30|] of the tree-level spinor-particle radiation rate for the special 
case of having only a single nonzero parameter a° > (i.e., case 2 of App. C of Ref. [l^ ) 
gives coefficient ^(qprim) = ('''/24) «° on the right-hand side of ([HJ. With ^ of order k > 0, 
the particle would slow down from an initial energy -Eprim ^ -f'thrcsh to the threshold energy 
i^thresh after having traveled a distance of the order of meters rather than parsecs (at least, 
for the values of k and -Eprim considered in this article). 

It may be worthwhile to repeat that expression ([8]) holds for an electric point charge, 
whereas we expect a form factor for a finite charge distribution with typical length scale 



a = he/ A] see also the discussion in Refs. 3l|, |3^. But this modification of the theory 



would primarily affect the total radiation rate and not the energy threshold ([9|) on which 
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our bounds are based, at least, for large enough cutoff A. Anyway, the theory considered in 
this article is precisely the one given in Sec. Ill Al without additional contact terms. 

Using the explicit (direction-dependent) result for the threshold energy ([9]), the con- 
dition -Eprim < -Ethresh cau then be written as the following upper bound on the deformation 
parameters of modified-Maxwell theory ([7]): 

i?(a° + a^ q^p,i^ + 5^'=q^p,i^q^J < [M^^im c'/ E^,,^)' , (10) 

with each index j and k summed over 1 to 3 and ramp function R{x) = {x + \x\ )/2. The 
parameters a^^ appearing in the argument of the ramp function on the left-hand side of (fTOj) 
are defined in ([5]), with 5^^ = —5^^ — 5^^ from the tracelessness condition ([3]). 

Observe that, as expected on general grounds, bound ffTU]) is invariant under a simulta- 
neous rotation of qprim and appropriate redefinition of parameters ■ ■ -a^. This becomes 
especially clear if the argument of the ramp function on the left-hand side of f IlOp is written 
as a° + (a ■ qprim) + (/3 ■ qprim) (7 ■ Qprim)' with two orthogonal 3-vectors j3 and 7 replacing 
the 3-tensor a^'' . In the same way, it is possible to absorb a parity-reflection of qprim by a 
change of sign of the three parameters entering ex., while leaving the other six parameters 
unchanged. 

In the following, we will use (fTOj) to obtain bounds on all nine deformation parameters 
«' from a sufficiently large number of UHECR events distributed over a large enough 
part of the sky (having primary energies En and directions q„, for n = 1, . . . , A^) and an 
appropriate value of Mprim- 



III. DATA SELECTION 



The Pierre Auger Collaboration has published a table of 29 UHECR events with energies 
above lOEeV = 10^^ eV 16|]. In addition, the atmospheric depth of the shower maximum. 



A^max, was determined for all of these events. Since the presently available Auger results 
indicate a mixed composition above 10 EeV js^l and our upper bound ffTOl) scales with M 



prim) 



we aim at selecting events initiated by light primaries primaries such as protons or helium 



nuclei. 
Ref. 



161: 



^or this reason, we apply the following selection criteria to the list of 29 events from 



(1) Xmax/(g cm-2) > 750 + 50 (log^o (E/EeV) - l); 

(2) > +2, where quantifies, in units of standard deviation, the difference between 



the observed X, 
of Ref. 



and the average X^^^,^ expected from primary photons [see Eq. (1) 



16| for the definition of A 



7 J' 



(3) primary energy between 10 and 30 EeV. 
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TABLE I: Simulated average shower- maximum atmospheric depth Xmax in units of g cm~^, with 
root-mean-square in brackets. Uncertainties are of the order of 3 gcm~^ for nuclear primaries and 
5 g cm^^ for photons. The calculations were performed with CORSIKA [sj] using two different 
hadronic interaction models, QGSJET 01 (abbreviated Q) and SIBYLL 2.1 (abbreviated 
S) [3]. The asterisk (*) indicates that, for primary photons above 10^^'^ eV, the Xmax values 
also depend on the direction of the event due to geomagnetic cascading (see, e.g., Ref. [37] and 
references therein). For instance, the X^^ax values at lO^'^ eV range between 940(85) and 1225(175) 
gcm~^, according to Ref. [ssl- 



logio(i^/eV) 


17.5 


18.0 


18.5 


19.0 


19.5 


20.0 


20.5 


21.0 


Q-P 


695(65) 


725(65) 


755(65) 


775(65) 


800(60) 


820(60) 


845(55) 


865(55) 


Q~He 


665(45) 


697(45) 


725(45) 


750(45) 


775(40) 


800(40) 


825(38) 


850(35) 


Q 






693(31) 


717(30) 


745(30) 


780(30) 


800(28) 




Q— Fe 


600(25) 


632(22) 


664(22) 


695(22) 


725(20) 


755(20) 


782(20) 


810(20) 


S-p 




740(65) 




800(60) 




860(55) 


885(50) 




S— Fe 




640(22) 




700(22) 




755(20) 


785(20) 




photon 






915(60) 


975(70) 


1075(95) 


(*) 


(*) 


(*) 



The first (and most selective) criterion aims at rejecting heavier nuclei. The chosen X^ax 
parametrization roughly follows the average Xmax values expected for helium nuclei. Thus, 
about half of the number of helium nuclei are rejected by this cut, while keeping the majority 
of protons. Oxygen and iron nuclei are rejected at a high level. In terms of the parameter 
Aprim, the selected events roughly have Aoxygcn < —1 and Ajron < —2, respectively. Table [T] 
gives some calculated values for Xmax as a function of primary energy and primary type 
(hydrogen, helium, oxygen, and iron nuclei and photon). 

The second criterion is to ensure no contamination from primary photons, as neutral 
photons would not emit vacuum Cherenkov radiation at tree level. In fact, all 29 events 
from the original table of Ref. 16|] already have A.^ > +2. 

The third criterion provides for a more homogeneous sample, but can certainly be relaxed 
in a future analysis. We have also accounted for the increased missing-energy correction in 
case of a nuclear primary, that is, the energies given in Ref. 16| (which refer to photon 
primaries) were increased by 7%, according to Ref. 39j. 

In total, 15 events remain which are listed in Table [Tll The corresponding event di- 
rections have not yet been published by the Pierre Auger Collaboration. Shown in Ta- 
ble [TTl are pseudo-random directions chosen from uniform distributions of right ascension 
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TABLE II: Selected Auger events from Ref. [Ig]: event identification number, primary energy E 
[EeV], shower-maximum atmospheric depth X^ax [gcm~-^], and pseudo-random event directions 
with right ascension RA' G [0°,360°] and declination 5' S [—70°, 25°]. The primes on the ID 
numbers are to emphasize the nonreality of these event directions, which can later be replaced by 
the measured values RA and 6 after their release by the Pierre Auger Collaboration (dropping the 
primes on the ID Nos.). 



ID No. 


E 


^max 


(RA' 


5') 


668949' 


18 


765 


(356, - 


-29) 


673409' 


13 


760 


(344, - 


-62) 


828057' 


14 


805 


(086, - 


-34) 


986990' 


16 


810 


(152,- 


-33) 


1109855' 


17 


819 


(280, - 


-30) 


1171225' 


16 


786 


(309, - 


-70) 


1175036' 


18 


780 


(228,- 


fl7) 


1421093' 


27 


831 


(079,- 


fl3) 


1535139' 


16 


768 


(006, - 


-62) 


1539432' 


13 


787 


(153,- 


-15) 


1671524' 


14 


806 


(028, - 


-63) 


1683620' 


21 


824 


(024, - 


-23) 


1687849' 


17 


780 


(031,- 


-23) 


2035613' 


12 


802 


(079, - 


-08) 


2036381' 


29 


782 


(158,- 


-03) 



RA' G [0°,360°] and declination 6' G [—70°, 25°], where the primes indicate the fictional na- 
ture of these directions. Note that the presently known UHECRs at energies E ^ 10 EeV are 
consistent with the hypothesis of large-scale isotropy; see, e.g., Refs. |l3L Il4j. l40l. l4ll. l42l. 1431. l44 1 
for a selection of the available data. 

All bounds of this paper are based on the values given in Table HIl but, for the reason 
given above, we expect that the same bounds are obtained from the actual event directions 
when they are made available by the Pierre Auger Collaboration. Just to be clear, the 
bounds of this article simply follow from having 15 UHECR events with more or less equal 
energies of order 10 EeV and more or less random directions over a significant part of the 
sky, the precise association of energy and direction being irrelevant. 
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The average energy of the sample given in Table HTl is {E) 11 EeV with an uncertainty 



of about 25% j45j. A larger uncertainty concerns the identity of the primary particle. We 
estimate an average mass for the selected events of (Mprim) ~ / x 5 GeV/c^, where the 
factor / G [0.2, 1.8] indicates the uncertainty (/ = 0.2 refers to a pure prot^ beam, / = 1 



to a mixed four- component primary composition which fits the Auger data [33|], and / = 1.8 
to the extreme case of no protons at all in the primary beam and a mixed three-component 
composition). In the next section, we will use an even larger value of /, namely / = 3.2, in 
order to obtain absolutely reliable bounds on the deformation parameters aK 



IV. UHECR CHERENKOV BOUNDS 

Following up on the remarks of the last paragraph of Sec. IIIBl it is now possible to 
determine numerically an exclusion domain in the a parameter space from the = 15 
events of Table [TTl and the Cherenkov threshold condition ffTU]) . This exclusion procedure 
only works for those parameters a for which the phase velocity of light is strictly less than 
the maximal attainable velocity c of the charged particle in theory ([7]), f phase < c. In fact, 
this phase- velocity condition corresponds to having a positive argument of the ramp function 
on the left-hand side of (fTOj) . Hence, the considered domain of parameter space is given by 

^causa/ = {5 G : V^eM3 (a° + P + > 0}, (11) 

with X an arbitrary unit vector in Euclidean 3-space, indices j and k summed over 1 to 3, 
and parameters a^^ defined by ([3]). As mentioned in Ref. [l^ and Sec. Ill A[ the condition 
"Wfront = "Wphase < c is, most likely, necessary for having a consistent and causal version of 
modified-Maxwell theory, hence the suffix "causal" in f|TT]) . The other suffix "(open)" in ffTTj) 
denotes the restriction of the causal domain by use of the open inequality "(a° + ■■■)> 0" 
instead of the closed inequality "(a'' + ■■■) > 0", because no Cherenkov limits can be 
obtained for fphase = c. 

A hypersphere of finite radius a in the subspace f lTTi) can be found numerically, so 
that for each a ^ the inequality (ITOl) is violated for at least one event from Table HTl 
The excluded domain of parameter space then corresponds to the region on or outside this 
hypersphere. Observe, namely, that the left-hand side of flTUl) . for a positive argument of 
the ramp function, increases by a factor A under the scaling of a ^ A a with A > 1, so that 
inequality (ITU]) is violated for A > 1 if it is already for A = 1. 

In this way, the experimental data of Table [III allow us to exclude the following region of 
parameter space at the two-cr level: 

Dt.,M,i = d'^ZI n {5 £ R» : |a| > a], (12a) 
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where |a| is the standard Euchdean norm and where the reference value for the mass 
of the primary charged particle has conservatively been taken equal to that of oxygen (as 
mentioned in Sec. IIIH most primaries in the sample are expected to be protons and helium 
nuclei). The 2a statistical error included in the number 3 x 10~^^ on the right-hand side of 
(]12bp is mainly due to the 25 % energy uncertainty of the UHECRs. 

The domain f|T2|) gives only part of the excluded region, because at certain points on the 
hypersphere the excluded region pushes in. An interesting question, then, is how far 
the excluded region extends inwards. A partial answer is given by the following result. It 
is possible to establish numerically the existence of a hyperball with finite radius b, so 
that for each a & the inequality (fTOj) holds for all events from Table [Tll Specifically, the 
following domain of parameter space is found to be allowed: 

^allowed = ^causal H {5 G : \d\ < b} , (I3a) 

6M1.5±0,7)xlO-(^^)\ (13b) 

with the approximate la error in fllSbp being mainly due to the energy uncertainty of the 



UHECRs and using the same reference value for the mass Mprim as in (]12b|) . The complete 
allowed region may very well extend beyond the radius b, but not beyond the radius a as 
given by fll2b|) . 



Returning to the exclusion domain (fT2|) . the corresponding 2a bound on |a| is given by 

where, as argued above, a value of 16 GeV/c^ for MpHm is entirely reasonable. With Lorentz 
invariance violated by the modified-Maxwell theory, it is important to specify the reference 
frame in which bound f|T^ holds, namely, the solar-system frame in which the cosmic-ray 
energies are measured. 

The new bound (1141) is consistent with earlier expectations, as given by the one-sided 



bound (C12) of Ref. [1^ in a somewhat implicit notation. Remark that our new bound on 
d is also "one-sided" (vacuum Cherenkov radiation occurs only if the light velocity is less 
than that of the charged particle), but, now, the one-sidedness is in a 9-dimensional space. 
In order to be as clear as possible, the exclusion domain f|T2l) has explicitly been given as an 
overlap with the domain (fTT!) defined at the beginning of this section and bound (IT^ has 
been specified to hold for a in the same domain. 

The domains f|T2l) and f|T3|) . together with the resulting upper bound (|T^ . have been 
calculated from the measured (unprimed) energies and pseudo-random (primed) directions 
shown in Table [ITl But, as mentioned in Sec. Illlt we expect that the same results are 
obtained if the measured (unprimed) event directions are used. 
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In closing, it may be helpful to rephrase our main result (HM as follows: considering 
deformation parameters a^, . . . with a corresponding phase velocity of light less than 
the maximal attainable velocity of charged particles, each of these nine parameters must 
separately have a modulus less than the value a given by (]12b|) at the 2a level. A more 
detailed analysis of these bounds is left to the future when more UHECR events become 
available. 



V. SUMMARY 

In this article, we have established an UHECR Cherenkov bound f[T^ on the nine de- 
formation parameters ([5]) of nonbirefringent modified-Maxwell theory ([I])-(l2]) coupled to a 
standard Dirac particle, with total action ([7j). Three remarks may be helpful to clarify the 
background and meaning of our result. 

First, the derivation of bound f|T^ relies, strictly speaking, on the assumption of large- 
scale isotropy of the arrival directions of 10-30 EeV cosmic rays, which is, however, supported 
by all the available data (see Sec. lIIII for further discussion). Second, this bound is essentially 
"one-sided" as it applies only to parameters belonging to domain (|TT1) . which may, however, 
correspond to a large part of the "physical domain" of the theory (see Sees. Ill Al and HVl for 
further discussion). Third, bound IHM holds for the theory as defined by the action ([7]) but 
also applies to a theory with additional Lorentz- violating parameters c^'^ in the Dirac sector 



2], provided the parameters k^'^ in bound f[T^ are replaced by the effective parameters 
K^j^ = kI''^ -2 c^^ 0, Q, 



Having clarified the background of our Cherenkov bound flMl) , it is to be remarked that the 
corresponding bound at the 10~^^ level for the spatially isotropic nonbirefringent deformation 
parameter a° improves significantly upon the direct laboratory bound at the 10^^ level 
[6(c)] or the indirect (electron anomalous-magnetic- moment) laboratory bound at the 10^® 
level [7(c)]. But also the other Cherenkov bounds at the 10~^^ level for the nonisotropic 
parameters a', for / = 1,...,8, improve considerably upon the direct laboratory bounds 
at the 10^^^ to 10^^^ levels jsl. While we note that a consistent comparison of bounds 
from laboratory and astrophysics data needs to be performed with great care (for example, 
Ref. [8(b)] allows for independent Lorentz- violating parameters in the photonic and fermionic 
sectors), the astrophysics bounds at the 10~^^ level are certainly indicative. 

With more and more UHECR events detected at the Pierre Auger Observatory over the 
coming years, these Cherenkov bounds at the 10~^^ level can be expected to improve on at 
least two grounds. First, the availability of a larger number of events and further observables 
in addition to X^ax will allow for more stringent cuts to obtain an essentially pure proton 
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sample, which would give a reduction factor of 16^ in (112bp .^ Second, higher energy events 
will become available for the analysis, perhaps up to values of order E' ~ 3 x 10^° eV 



(the current highest known energy 46|]), which would give a further reduction factor of 
approximately 10'^ in (]12b|) .^ Combined, the Cherenkov bounds on the nine deformation 
parameters a' can perhaps be reduced by a factor of order 10^ to the 10~^^ level in the next 
decade. 
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NOTE ADDED IN PROOF 

The Pierre Auger Collaboration has recently published a list of 27 UHECR events with 
energies above 57 EeV and accurate event directions, which were recorded over the period 



1 January 2004 to 31 August 2007 49|]. In order to provide better coverage of the northern 



celestial hemisphere, we add a 210 EeV event from the AGASA array [50i] and the previously 



mentioned 320 EeV event from the Fly's Eye detector [46|], both with reasonably accurate 
arrival directions. The observed energies and directions of these 29 events can now be used 
to sharpen bound f|T^ in terms of the square of the number a given by (]12bp . 

In the absence of specific information on the chemical composition of all 29 events, we 
simply take Mprim ^ 16 GeV/c^. From these particular 29 UHECR events with an estimated 
energy uncertainty of the order of 25 %, the improved value for a is found to be given by 

A ^ a.,„ ^ 2 X 10-1^ ( ' ' (15) 



16 GeV/c2 

which gives the 2a excluded region -D^^j^^g^ as defined by I12a[ With the 29 event directions 
known explicitly, there is no assumption entering the corresponding 2a bound f[T^ with 
on the right-hand side, apart from the uncertainty in the primary mass and the caveat 
mentioned in Sec. IIIBI 



"'^ If a sufficient number of 10 EeV cosmic-ray electrons could somehow be identified, there could be a further 

reduction factor in (|12bp of approximately (938/0.511)^ w 3 x 10® compared to the proton case. 
^ The steeply falling flux spectrum and the expected suppression of protons from cosmologically distant 



sources with energies above 6 x 10^^ eV 47|, |48| present an experimental challenge for pushing the energy 
to substantially higher values. 
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